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Motivation 


- Perform detailed time-accurate scale-resolving simulations of practical, 
complex, compressible flows 


¢ High-Reynolds-number separated flows involving large-scale 
unsteadiness, where RANS models are unreliable 
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Approach 


- Turbulent flows involve a large range of spatial and temporal 
scales which need to be resolved 


¢ Efficient algorithms and implementation necessary for wall- 
resolved high Reynolds number flows 


¢ Numerical method must be capable of handling complex 
geometry 


¢ Numerical method must be “‘robust”’ 


Developed higher-order space-time discontinuous-Galerkin 
spectral element framework 


o 


Approach 


¢ Gradient computation needed for error-estimation, adaptation, 
design, sensitivity analysis, etc. 


- Tangent and adjoint methods have been successfully applied to 
a variety of steady and unsteady flows 


- High-fidelity simulations we are targeting are chaotic 
-Can traditional tangent/adjoint methods work? 


¢Develop efficient implementation of tangent and adjoint 
method in space-time discontinuous solver 


¢ Assess the applicability of traditional adjoint and tangent 
methods for chaotic flows 


L. Diosady 03-01-17 ehkhy 4 


Outline 


¢ Space-time DG formulation 
- Discrete primal formulation 
- Discrete tangent and adjoint formulations 
* NACAOOI2 
* Flow sensitivity with increasing Reynolds number 
- Properties of adjoint for chaotic flows 
-T106a LPT 
¢ Adjoint solutions corresponding to a practical simulation 
- Summary/Outlook 
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Space-Time Discontinuous-Galerkin (DG) asa 


Formulation 


Compressible Navier-Stokes Equations: = 


Ot 
¢ Space-Time Discontinuous-Galerkin Discretization: 


+V-F(u, Vu) = 0 


- Entropy variables: Hughes (1986) 


SPD Sym  SPSD a1 
= 


Aov 4 ++ Ayu j a (Ki 50,5) = (0) -— 
¢ Inviscid Flux: Ismail & Roe (2009) 


¢ Viscous Flux: Interior penalty method with penalty parameter 
given by 2nd method of Bassi and Rebay (2007) 


- Integrals evaluated using numerical quadrature with 2N points 

¢ Discrete entropy stability: Barth (1995) 
(ps) ,¢ + 7 
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Discrete Formulation S 
- Discrete system solved each time-slab: 


R”(u™,w") + G"(u, w") =0 


¢ where: 


ll -S{ff- (Siu + Vu) tf f or ne [we ue)! 
G"(u"-1,w”) = + | -w(Eu(t") f 


RK 


- beginning/end of time-slab 


w(t) = I,w™ w(t?tt) = Iw” 
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Discrete Tangent formulation 


° Output: J(u; cu) 


* where © 1S a parameter (1.e. angle of attack, Reynolds number etc.) 


Compute sensitivity of output to parameters: 


dJ_aJ  aJdu _ as , al aR 
da Oa Ouda Oa OROAa 
- langent equation: 
OR” n n OG” n—1 N\ OR” 
5,7 (OU ,Ww ) + 5-a=7 (Ou ,W i = 
¢ Matrix form: 
OR" 0 Sun! Onn 
aG” aR” n OR” 
0 OU 
Oun—1t Ou” Ou” 
0 AG't! aRrt! oyrtl ORT! 
Ou” Our Our 
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Discrete Adjoint Formulation 


- Lagrangian: L(u, wb; a) = J(u;a) +" R(u; a) 


* Stationarity of Lagrangian: 


as* OR Sh id of 
s\£L= | — . OU \0a=0 
2 & ? Ou ) (S| +e ‘| ° 
Aqyjoint ne 
¢ Adjoint equation: 
OR” n n OG” n—1 N\ OJ” n 
Ou” (w mI ) + Our-1 (w ,w ) — Ayn (ew ) 
¢ Matrix Form: 
O n—-1d n LT n—l 
en cee 0 ue —— 
_ peti 7° er ones 
0 0 dynrtl Our? 
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Primal Solver: Implementation Details 


¢ Efficient implementation of higher order DG 
- Tensor-product basis 
- Take advantage of hardware (SIMD/optimized kernels) 
¢ Jacobian-free Approximate Newton-Krylov solver 
- Tensor-product based ADI-preconditioner 


¢ Primal Residual Evaluation: 


Optimized 


1. Evaluate state/gradient at quadrature points aaa 
sum-factorization 


2. Evaluate flux at quadrature points 


3. Weight fluxes with gradient of test functions Wertereed Warnick 
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Tangent/Adjoint: Implementation Details © 


¢ Reuse optimization from primal for tangent and adjoint 


- Tangent Residual Evaluation: 


1. Evaluate state/gradient and update/gradient at quadrature points 


2. Evaluate linearized flux at quadrature points — Optimized 


slimaacior ization 


3. Weight fluxes with gradient of test functions 
¢ Adjoint Residual Evaluation: 

4. Evaluate state/gradient and adjoint/gradient at qué &points 

5. Evaluate adjoint flux at quadrature points 


Vectorized Kernels 
6. Weight fluxes with gradient of test functions 
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NACA0012, «= 10 Nasa] 


- At high angle of attack, flow over NACAOO12 airfoil exhibits 
vortex shedding which become chaotic with increasing Reynolds 
number. (Pulliam 1993) 


¢ Examine primal and adjoint solutions 
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NACA0012, Re = 800, a = 10 Nasa] 


Adjoint Magnitude 
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- Unsteady shedding gives periodic output signal 


¢ Adjoint solution also periodic (after initial transient) 
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NACAOO12, Re = 1600, a = 10 
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Adyjoint Magnitude 


¢ With increasing Reynolds number force has multiple frequencies 


* Adjoint solution still essentially appears periodic 
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NACA0012, Re =2400, a = 10 Nasa] 
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¢ With increasing Reynolds number simulated flow become chaotic 


- Adjoint solution begins to grow unboundedly 
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NACAO0012, Re =800-2400,a=10 @& 


Adjoint Magnitude 


-Time 


Adyjoint Magnitude 


¢ With increasing Reynolds number simulated flow become chaotic 


- Adjoint solution begins to grow unboundedly 
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NACA0012, Re =800-2400,a=10 @& 
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Adyjoint Magnitude 
- Solution is in fact chaotic at Re = 1600, but growth rate is much 
slower than at Re = 2400 


¢ Windowing approaches may be successful at Re = 800, 1600 
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T106c Low Pressure Turbine Nasa] 


-Re = 80,000, Minflow = 9.243, alpha = 32.7, M..:, = 0.65 
¢ Periodic BCs in span-wise and pitch-wise directions 
- No free-stream turbulence 


- Spanwise domain is 20% of chord 
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Primal Solution asa 


* Nearly steady flow 
upstream and over first 
2/3 of blade 


- Separation leading to 
transition/vortex 
shedding on suction- 
side of blade 


- Fully turbulent wake 
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Sensitivity to Inflow Boundary Condition Nasa] 


* Modify inlet flow angle from a = 32.7 to a = 32.701 


Convective disturbance Domain flow- 
hits leading edge through time 
1.5 | 3 
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Sensitivity to Inflow Boundary Condition Nasa] 


Convective disturbance 


hits leading edge 
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Adjoint of mean Axial Force 


- Output is integrated axial force 


— 4. rt 
f= 7 | F., (u(r) )dr 


- Also define output without 
temporal normalization 


Range: [-le6, |e6] 


L. Diosady 03-01-17 ehkhy vy) 


Sensitivity computed using adjoint oS 
AJ(t) = J(t;a+ Aa) — J(t;a) = / F,,(u(t;a + Aa)) — F, (u(t; a))dr 


~ | W(7;t,a)* R(u(r);a + Aa) 


lJ(at+t A a)-J(a)! 
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|[Axial Force Differencel 
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Sensitivity computed using adjoint eS 


AF,,.(7) = Fy (u(t;a+ Aa)) — F,(u(t; a)) 
2? & W(t —7;t,a)’ R(u;a+ Aa) 
¢ Approximation holds 


since flow upstream of 
blade 1s essentially time- 


independent 
— _ = - Adjoint correctly captures 
— join —(). a . 
Adjoint T=0 2 sensitivity in part of flow 
- = = Adjoint T=0.3 upstream of separation 


Adjoint T=0.4 
Adjoint T=0.7 


Time 
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Sensitivity computed using adjoint ] 


Time 


- Sensitivity computed 


using adjoint only valid 
for very short time 
windows 


- Adjoint computed using 


long time window blows 
up 


- Sensitivity computed 


using short time window, 
not representative long 
time behaviour 


ehkhy 
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Graphical demonstration of concepts 


¢ Adjoint correctly captures 
sensitivity in part of flow 
upstream of separation 


- Sensitivity computed using 
adjoint only valid for very 
short time windows 


- Sensitivity computed using 
short time window, not 
representative long time 
behaviour 
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Adyjoint-based Error Estimation S 


¢ Estimate error using dual-weighted 
residual method (Becker & 
Rannacher 1995) 


e= J(u) —J(un) © Ra(ug, wn) 


Primal Solution 


¢ Localize error 
€;, = Ry(ug,Vnalx) 


- Flag elements with largest error for 
refinement 


Element-based error-indicator 
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Adyjoint-based error indicator Nasa] 


log(error) 


¢- Unbounded adjoint not useful for error 
estimation 


- Estimate is orders of magnitude larger than 
actual signal 


- Error localization simply flags regions where 
adjoint is large 
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Adjoint growth with mesh resolution 


- Refined mesh has essentially 
double mesh resolution near 
separation region 


¢ Increase mesh resolution 
results in faster growth of 
adjoint (1.e. larger Lyapunov 
exponent) 


- Adaptation mechanism is not 2 
convergent 


Time 
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Summary Nasa] 


- Presented space-time adjoint solver for turbulent compressible 
flows 


¢ Confirmed failure of traditional sensitivity methods for chaotic 
flows 


- Assessed rate of exponential growth of adjoint for practical 3D 
turbulent simulation 


- Demonstrated failure of short-window sensitivity 
approximations. 
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Questions?’ 


Outlook/Future Work: 


- Lyapunov exponents, least-square shadowing and beyond... 
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